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Abstract 



We deal with a nonconvex and nonlocal variational problem coming from thin-film micro- 
magnetics. It consists in a free-energy functional depending on two small parameters e and rj 
and defined over vector fields m : Q C R 2 — > S 2 that are tangent at the boundary dSl. We 
are interested in the behavior of minimizers as £,77 — !> 0. They tend to be in-plane away from 
a region of length scale e (generically, an interior vortex ball or two boundary vortex balls) 
and of vanishing divergence, so that S 1 — transition layers of length scale n (Neel walls) are 
enforced by the boundary condition. We first prove an upper bound for the minimal energy 
that corresponds to the cost of a vortex and the configuration of Neel walls associated to 
the viscosity solution, so-called Landau state. Our main result concerns the compactness of 
vector fields {m e ,^} £ ^4o of energies close to the Landau state in the regime where a vortex 
is energetically more expensive than a Neel wall. Our method uses techniques developed for 
the Ginzburg-Landau type problems for the concentration of energy on vortex balls, together 
with an approximation argument of S 2 — vector fields by S 1 — vector fields away from the vortex 
balls. 

AMS classification: Primary: 49S05, Secondary: 82D40, 35A15, 35B25. 
Keywords: compactness, singular perturbation, vortex, Neel wall, micromagnetics. 

1 Introduction 

In this paper, we investigate a common pattern of the magnetization in thin ferromagnetic films, 
called Landau state, that corresponds to the global minimizer of the micromagnetic energy in a 
certain regime. For that, we focus on a toy problem rather than on the full physical model: 

Let n C R 2 be a bounded simply-connected domain with aC 1,1 boundary corresponding to the 
horizontal section of a ferromagnetic cylinder of small thickness. Due to the thin film geometry, 
the variations of the magnetization in the thickness direction are strongly penalized. It motivates 
us to consider magnetizations that are invariant in the out-of-plane variable, i.e., 
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and they are tangent to the boundary dfl, i.e., 



m'-v = Q on <9ft, (1) 

where m! — (mi, m%) is the in-plane component of the magnetization and v is the normal outer 
unit vector to <9f2. We consider the following micromagnetic energy functional: 

E e Jm)= [ \\7m\ 2 dx + \ [ mldx+- [ \\V\~ 1/2 (V ■ m')\ 2 dx, 
Ju £ Jn V Jut 

where e and rj are two small positive parameters (standing for the size of the vortex core and the 
Neel wall core, respectively). Here, x = (x\,X2) are the in-plane variables with the differential 
operator 

V = (d xl ,d X2 ). 

The first term of E £ ^(m) stands for the exchange energy. The second term corresponds to 
the stray-field energy penalizing the top and bottom surface charges 7713 of the magnetic cylinder, 
while the last term counts the stray-field energy penalizing the volume charges V • m' where we 
will always think of 

to' = to'Iq 

as being extended by outside fl. For more physical details, we refer to Section [3l 

Note that the non-local term in the energy is given by the homogeneous if -1 / 2 — seminorm of 
the in-plane divergence V • ml that writes in the Fourier space as: 

2 

||V-m'|||_ 1/2(R2) = / iVr^V-m') dx:= [ ±.\F(y . m ')\* (2 ) 

Also observe that the boundary condition ([1} is necessary so that @ is finite since 

V • m' = (V • to')1q + (to' • v)l an in R 2 

(see Proposition [5] in Appendix). 

We are interested in the asymptotic behavior of minimizcrs of the energy E e ^ in the regime 

e«l and i] <C 1. 

The main features of this variational model resides in the nonconvex constraint on the magnetiza- 
tion \m\ = 1 and the nonlocality of the stray- field interaction. The competition of these effects with 
the quantum mechanical exchange effect leads to a rich pattern formation for the stable states of 
the magnetization. Generically, a pattern of a stable state consists in large uniformly magnetized 
regions (magnetic domains) separated by narrow smooth transition layers (wall domains) where the 
magnetization varies rapidly. The characteristic wall domains observed in thin ferromagnetic films 
are the Neel walls (corresponding to a one-dimensional in-plane rotation connecting two directions 
of the magnetization) together with topological defects standing for interior vortices (called Bloch 
lines) and micromagnetic boundary vortices. 

The existence of line singularities at the mesoscopic level of the magnetization in thin films can 
be explained by the principle of pole avoidance. For this discussion, we first neglect the exchange 
term in E e , r/ . The stray-field energy will try to enforce in-plane configurations, i.e., to 3 = in f2, 
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together with the divergence- free condition for m', i.e., V • mf — in ft. Together with (p}, we 
arrive at 

|m'| = 1, V • mf = in O and mf ■ v = on 90. (3) 

We notice that the conditions in ([3]) are too rigid for smooth magnetization mf . This can be seen 
by writing m! = V ip with the help of a "stream function" ip. Then up to an additive constant, 
([3]) implies that ip is a solution of the Dirichlet problem for the eikonal equation: 

\Vtp\ =1 in Q and V = on dSl. (4) 

The method of characteristics yields the nonexistence of smooth solutions of But there are 
many continuous solutions that satisfy (j4j away from a set of vanishing Lebesgue measure. One 
of them is the "viscosity solution" given by the distance function 

tp(x) = dist (x,dCl') 

that corresponds to the so-called Landau state for the magnetization m! . Hence, the boundary 
conditions (fTJ) are expected to induce line-singularities for solutions mf that are an idealization 
of wall domains at the mesoscopic level. At the microscopic level, they are replaced by smooth 
transition layers, the Neel walls, where the magnetization varies very quickly on a small length 
scale 1]. Note that the normal component of mf does not jump across these discontinuity lines 
(because of (J3j) ) ; therefore, the normal vector of the mesoscopic wall is determined by the angle 
between the mesoscopic levels of the magnetization in the adjacent domains (called angle wall). 
Now, taking into account the contribution of the exchange effect, the energy scaling per unit length 
of a Neel wall of angle 29 (with 9 € (0, |])) is given in DeSimone, Kohn, Miiller& Otto [f\, Ignat 
&Otto pi] (see also Ignat [8]): 

tt(1 -cos#) 2 + o(l) 
?7|log7?| 

The formation of interior or boundary vortices is explained by the competition between the 
exchange energy and the penalization of the 7713— component for configurations tangent at the 
boundary. Indeed, there is no S 1 — configuration that is of finite exchange energy and satisfies (fl}. 
There are only two possible situations: If m! does not vanish on dQ, than (TIJ implies that m! 
carries a nonzero topological degree, deg(m',<9f2) — ±1. In this case, we expect the nucleation 
of an interior vortex of core-scale e. The scaling of the vortex energy is related to the minimal 
Ginzburg-Landau (GL) energy (see Bethuel, Brezis & Helein [T]): 

min / q £ (m') dx = (2tt + o(l))\ logsl as e — > 0, (6) 

m'—u 1 - on dQ 

where the GL density energy is given in the following: 

5£ (m')HVmf + l(l-|mf) 2 . (7) 

(Here, we denote v = (— v?,,v\)f) The second situation consists in having zeros of mf on the 
boundary. Therefore, we expect that boundary vortices do appear. Roughly speaking, they corre- 
spond to "half of an interior vortex where the vector field m! is tangent at the boundary; therefore 
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they arc different from the micromagnetic boundary vortices analyzed by Kurzke[14] and Moser 
[IB] (see details in Section [3]) . Remark the importance of the regularity of dil in estimate ^ . In 
fact, if dfl has a corner and the boundary condition m! = v on dVL in ([6]) is relaxed to (JTJ) , then 
estimate (O does not hold anymore, it depends on the angle of the corner (see Proposition Q] and 
Remark [5]). Therefore, at the microscopic level, topological point defects do appear in the Landau 
state pattern and are induced by ((TJ) . 

The aim of the paper is to show compactness of magnetizations of energy E EtV close to the 
Landau state in order to rigorously justify the limit behavior ([3]): the delicate issue consists in 
having the constraint \m\ = 1 conserved in the limit. For that, we have to evaluate the energetic 
cost of the Landau state. We expect that the leading order energy of a Landau state is given by 
the topological point defects and Neel walls. The Landau state configuration consists in several 
Neel walls and either one interior Bloch line or two "half Bloch lines placed at the boundary of 
the sample J7. Therefore, by ([5]) and (j6]), we expect that the energy of the Landau state has the 
following order: 



2tt| loge| 



.4 



7/|log7/|' 

for some positive A > depending on the length and angle of Neel walls. 



(8) 



2 Main results 

First of all, we want to rigorously prove the upper bound ((SJ for the Landau state. Our result 
gives the exact leading order energy of the Landau state in the case of a domain f2 of a " stadium" 
shape (see Figure [T|) . Note that the Landau state of a stadium consists in a single Neel wall of 




Figure 1: Stadium 

180° (in our example, the length of the wall is equal to 2, so that A = 2ir in |8])). 
Theorem 1 Let Q — £!i U SI? U 1^3 be the following "stadium" shape domain: 



fii = {x = (xi,x 2 ) £ I 

n 2 = (-i,i)x (-1,1) 
tt 3 = {x = (xi,x 2 ) e I 



-(1,0)| < l.ii > 1}, 



(-1,0)1 < < -1}. 



In the regime e <C r\ -C 1, there exists a C 1 vector field m e . xi : ft — > S 2 that satisfies JTJ) and 
E e , v (m e<v ) < 27r|loge| + 2?r + °W as n | q. 

7/| log 7^1 



(9) 



Observe that the vortex energy in the above estimate is relevant only if a vortex costs at least 
as much as a Neel wall, i.e., ^ < | loge| (otherwise, the vortex energy would be absorbed by 
the term o( > ^ , j ). This regimes leads to a size e of the vortex core exponentially smaller than 
the size of the Neel wall core 77 (see Remark [T]). 

Notation: We always denote a <C b if ^ —> 0; also, a < b if a < C& for some universal constant 
C > 0. 

Now we state our main result on the compactness of the S 2 — valued magnetizations that have 
energies near the Landau state. The issue consists in rigorously justifying that the constraint 
\m\ = 1 is conserved by the limit configurations as e, n — > 0. The regime where we prove our result 
corresponds to the case where a topological defect is energetically more expensive than the Neel 
wall, that is coherent with the regime where © holds. 

Theorem 2 Let a G (0, ^) be an arbitrary constant. We consider the following regime between 
the small parameters e, r\ <C 1: 

e 1/2 < V, (10) 
log|loge|<^^. (11) 
For each e and rj, we consider C 1 vector fields m £ ^ v : — > S 2 that satisfy |T]) and 

f <27ra|loge| (12) 

Ee, v (m e , n ) -27r|loge| 1 1 

< -n — r (13) 

r)\\ogr)\ 

Then the family {m e . r,} £ ^^o *s relatively compact in L 1 ^^ 2 ) and any accumulation point m : 
SI — > S 2 satisfies 

TO3 = 0, \m\ = 1 a.e. in Q and V • m! = distributionally in M 2 . (14) 



The proof of compactness is based on an argument of approximating S 2 — valued vector fields 
by 5' 1 -valued vector fields away from a small defect region. This small region consists in either one 
interior vortex or two boundary vortices. The detection of this region is done in Theorem [3] and 
uses some topological methods due to Jerrard [12] and Sandier [18] for the concentration of the 
Ginzburg-Landau energy around vortices (see also Lin [T5], Sandier & Serfaty [TS] )■ Away from 
this small region, the energy level only allows for line singularities. Therefore, the compactness 
result for S 1 — valued vector fields in pj] applies. 

Let us discuss the assumptions (fT0]l . ([TT]). (1P2]) & (fT3|) . Inequality (IT3]) assures that cutting 
out the topological defect (one vortex or two boundary vortices) , the remaining energy rescaled at 
the energetic level of Neel walls is uniformly bounded. Inequality (IT2"]) together with the choice of 
a < I mean that the energy cannot support three " half interior vortices and is precisely explained 
in Theorem [3] below. Inequality (|11[) is imposed due to our method to detect a boundary vortex: it 
leads to a loss of energy of order 0(log | log s\) with respect to the expected half energy of a interior 
vortex, i.e., tt\ log e\ (see Theorem[3]and Proposition!!]) . This amount of energy could leave room for 
configurations of Neel walls that may distroy the compactness of |m'| = 1. Therefore, to avoid this 
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scenario, is imposed. The regime (1101) is rather technical: it is needed in the approximation 
argument of S 2 — valued vector fields by S 1 — valued vector fields away from the vortex balls. In fact, 
starting from the values of m! on a square grid of size , the approximation argument requires 
zero degree of m' on each cell, leading to the condition ft < 1 — a (see Lemma [2]); furthermore, the 
condition e' 3 < rj is needed in order that the approximating S 1 — valued vector fields induce a stray 
field energy of the same order of m! (see (|77| ). Therefore, (fTt)|) can be improved to a larger regime 

£ ^ ?5 V f° r anv /3 < l — a 

as presented in the proof (Theorem [2] is stated for the value (3 — 1/2 which is the universal choice 
for every a < 1/2). However, this slightly improved condition is weaker than the complete regime 
implied by (fT2|) as explained in the following remark. 

Remark 1 Any limit configuration m' satisfies (fT4|) . If is a bounded simply-connected domain 
different than discs, m! has at least one ridge (line-singularity) that corresponds to a Neel wall. 
Therefore, the minimal energy verifies minjjj E e-ri — 2ir\ \oge\ > , 1 * g , . Combining with (|12[) . it 
follows that 

-n — i ~ l lo s £ l; 

l 

in particular, e < e 51 lo sii , i.e., the core of the vortex is exponentially smaller than the core of the 
Neel wall. However, in the proof of Theorem [21 this much stronger constraint with respect to (fT0|) 
is not needed. 

We prove the following result of the concentration of Ginzburg-Landau energy around one 
interior vortex or two boundary vortices for vector fields tangent at the boundary: 

Theorem 3 Let a € (0, 5) and O C M 2 be a bounded simply- connected domain with a C ' bound- 
ary. There exists Eq — £q((X, dQ,) > such that for every < e < Eq, if ml : — > B 2 is a C 1 vector 
field that satisfies ([I]) and 

/ g £ (m')dx< 27r(l + a)|log£|, (15) 
Jn 

then there exists either a ball B{x\,r*) C £1 (called vortex ball) with r* = | lo g £ | 3 and 

I g e {m')dx>2K\\og-\-C, (16) 

JB{xl,r*) £ 

or two balls B(x2,r*) and B{x't > ,r*) (called boundary vortex balls) with X2,x% £ dfl and 

I g e (m')dx>27r\\og-\-C, (17) 

J(B(x*,r*)UB(xl,r*))nn £ 

where C = C{a,d£l) > is a constant depending only on a and on the geometry of dQ. 

The condition a < 1/2 is needed in our proof. In fact, if no topological defect exists in the 
interior (in which case, condition (fTJ induces boundary vortices), we perform a mirror- reflection 
extension of m! outside the domain. Roughly speaking, the GL energy in the extended domain 
doubles, i.e., it is of order 27r(2 + 2a)| loge| and the degree at the new boundary is equal to two; 
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in order to avoid the formation of three interior vortices in the extended region, we should impose 
2 + 2a < 3, i.e., a < 1/2. 

Notice that the Ginzburg-Landau energy concentration for a boundary vortex in (fl7|l has a 
cost of order ir\ loge\ — Clog | loge| provided that the boundary has regularity C 1 * 1 . We conjecture 
that the same energetic cost for a boundary vortex holds true if the boundary has regularity C 1 '^, 
/3 G (0, 1). However, if the boundary regularity is only C , then the energetic cost of a boundary 
vortex may decrease to (ir— log | Y og£ | )l l°g e l where C > is a universal constant. This indicates that 
the loss of energy of order log | loge| in (|17[) could occur for boundary vortices for C 1 '^ boundary 
regularity and the order of this loss increases to i g°iog £ | for C 1 boundaries as f3 — > 0. This claim 
is supported by the following example for a C 1 boundary domain: 

Proposition 1 We consider in polar coordinates the following C 1 domain f2 = {(r, 0) : r G 
(0, t^j), |6>| < 7(7*) = f — log log } ' F° r ever y < e < 1, there exists a C 1 — function m' e : 
f2 n Bi/200 — > R 2 that satisfies |T]) on d£l n B1/200 

/ g e (m' e )dx< (tt- - — ,)|loge|, 

JnnB 1/200 log I log e| 

where C > is some universal positive constant (independent of e). 

The outline of the paper is as follows. In Section [3J we present the physical context of our toy 
problem. In the next section, we recall two results that we need for the proof of our results: a 
compactness result for S 1 — valued magnetizations and the concentration of the Ginzburg-Landau 
energy on vortex balls. In Section [5j we prove Theorem [3] and Proposition [TJ In Section [6l we 
give the proof of our main result in Theorem [2j In Section [3 we show the upper bound for the 
stadium domain stated in Theorem [1] In Appendix, we prove that (fTJ is a necessary condition for 
our configurations to have a finite stray field energy. 



3 Physical context 

In this section we explain the physical context of our model in thin-film micromagnetics. We 
consider a ferromagnetic sample of cylinder shape, i.e. 

uj = J x (0,t) 

where u/ C R 2 is the section of the magnetic sample of length i and t is the thickness of the 
cylinder. The microscopic behavior of the magnetic body is described by a three-dimensional unit- 
length vector field m = (to', 7713) : u> — > S 2 , called magnetization. The observed ground state of 
the magnetization is a minimizer of the micromagnetic energy that we write here in the absence 
of anisotropy and external magnetic field: 

E 3d (m) = d 2 I |(V '^H 2 dxdz + f \{V,^-)U{m)\ 2 dxdz. (18) 

Ju UZ J R 3 oz 

The parameter d of the material is called exchange length and is of order of nanometers. The 
stray-field potential [/(to) : M 3 — > M. is defined by static Maxwell's equation in the weak sense: 

/ (V, ^-)U(m) ■ (V, |-)C dxdz = f (V, • (ml u ) C dxdz, for every £ G C^°(R 3 ). (19) 
J R 3 Oz dz J R3 dz 
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Instead of the three length scales £, t and d of the physical model, we introduce two dimensionless 
parameters: 

d , d 2 

e := - and r, := - 

(standing for the size of the core of the Bloch line and the Neel wall, respectively). 

Thin-film reduction. Wc consider the thin-film approximation of the full energy (fT8|) in the 
following regime: 

e < v < 1 (20) 

(equivalently, t <C d <C £)■ The assumption t <C d implies that in-plane transitions (Neel walls) are 
preferred to out-of-plane transitions (asymmetric Bloch walls) between two mesoscopic directions 
of the magnetization (see Otto [H]). The hypothesis d <C I assures that constant configurations 
in general are not global minimizers (see DcSimone [4]). 

The main issue is the asymptotic behavior of the energy in the regime of thin films. We first 
nondimensionalize in length with respect to £, i.e. (x, z) — (f,f), O = t-, fh(x,z) = m(x, z), 
U(rh)(x, z) = jU(m)(x, z) and then we renormalize the energy E 3d (fh) = ^E 3d (m). Omitting 
the " , we get 

E zd {m) = \ I |(V,^-)m| 2 dxdz+\ [ |(V,|-)[/(m)| 2 dxdz. (21) 

e 2 Jn X (o,4) dz e 4 J R3 dz 

In the regime (f20|). the penalization of exchange energy enforces the following constraints for the 
minimizers: 

2 

(a) m varies on length scales ^> ~. 

(b) m = m{x), i.e. m is z— invariant. 

With these assumptions, (f2Tj) can be approximated by the following reduced energy E red (see 
DeSimone, Kohn, Miiller & Otto 6 ], Kohn & Slastikov [13]): 

E red (m) = [ |Vm| 2 dx 
Jn 



- 2 L mldX ^^^ lj m '' V)2ml+ ^•^lll-W (22) 



The above formula follows by solving the stray field equation (jl9)) in the regime ([20]) : indeed, for 
z— invariant configurations m, the Fourier transform in the in-plane variables x = (xi,X2) turns 
(fT9|) into a second order ODE in the z— variable that can be solved explicitly (see [13], [9]). Then, 
due to the above assumption a) and to the regime (|20|). the stray-field energy asymptotically 
decomposes into three terms as written in (|22|) : the first term in (|22|) is penalizing the surface 
charges 7713 on the top and bottom of the cylinder, a second term counts the lateral charges m! ■ v 
in the L 2 — norm, as well as the third term that penalizes the volume charges (V-to') qc :— (V-77i')lf2 
as a homogeneous H^ 1 / 2 — seminorm. In fact, the last term corresponds to the stray-field energy 
created by a three-dimensional vector field ft. a c(w) defined as 

hac(m) = (V, ^-)U ac (m) : R 3 -> R 3 , 
dz 
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that satisfies: 



(V, — )U ac {m) ■ (V, ^-)C dxdz 



(V • m') ac C ^, for all C € C C °°(R 3 ). 



Then one has 



|ft-ac(m)| 2 <ia;dz = -||(V • m) ac \\ 2 ^_ 1 



/ 2 ( 



(23) 



Note that if JT]) holds (i.e., no lateral surface charges), then (V • m) ac = V • (mljj) and therefore, 
h a c{m) induces the stray field energy given by ©. In fact, corresponds to the minimal stray 
field energy in thin films. More precisely, a stray field h — {h! ,h$) = (hi, li2, h^) : R 3 — > R 3 is 
related to the magnetization m : Q — > S 2 via the following variational formulation: 

/ (ti -VC + hs^f) dxdz = [ (V-m'dx, VCeC c °°(R 3 ), (24) 
Jr 2 xr V vzj J R 2 

(As before, ml = m'lji and m satisfies 



where z denotes the out-of-plane variable in the space 
([!]).) Classically, this is, 



_ [h 3 ] = -V • m! 







m 

on 



\ (M 2 x {0}), 
2 x {0}, 



where [/13] denotes the jump of the out-of-plane component of h across the horizontal plane R 2 x {0}. 



Then §Z§ can be expressed as: 

/ |Vr 1/2 (V-m') 



dx = 2 min / \h\ dxdz 

7iwith|24j Jr2 xR 



Therefore, h ac (m) is a minimizing stray-field (of vanishing curl) associated with the stray field 
potential U ac (m). 

In our regime (|20|) . there are three different structures that typically appear: Neel walls, Bloch 
lines and micromagnetic boundary vortices. We explain these structures in the following and 
compare their respective energies. As we already mentioned, a fourth structure, the asymmetric 
Bloch wall, can appear in thicker films but we do not discuss it here since the asymmetric Bloch 
wall is more expensive than a Neel wall if t <C d. 

Neel walls. The Neel wall is a dominant transition layer in thin ferromagnetic films. It is 
characterized by a one-dimensional in-plane rotation connecting two (opposite) directions of the 
magnetization. It has two length scales: a small core with fast varying rotation and two logarith- 
mically decaying tails. In order for the Neel wall to exist, the tails are to be contained and we 
consider here the confining mechanism of the steric interaction with the sample edges. Typically, 
one may consider wall transitions of the form: 



m = (mi, 7712) : M — > S 1 and m(±£) 



cosw 
, ±sin6» / 

with 9 £ [0, -|) (see Figure 0, whereas the reduced energy functional is: 

dm 



for ±t> 1, 



E red (m) 



dxi 



' dxi 



2r] 







d 


1/2 


f 




dxi 


mi 



dx\. 
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Figure 2: Neel wall of angle 26 confined in [—1, 1]. 

As -q — > 0, the scale of the Neel core is given by \xi\ < w cor e — 0(rj) while the two logarithmic 
decaying tails scale as w core < < w t aii — 0(1). The energetic cost (by unit length) of a Neel 
wall is given by 

£i red (Neel wall) = 0( , - 1 , ) 
?7| log?7] 

with the exact prefactor ir(l — cos#) 2 /2 where 26 is the wall angle (see e.g. [5]). 

Bloch line. A Bloch line is a regularization of a vortex on the microscopic level of the magneti- 
zation that becomes out-of-plane at the center. The prototype of a Bloch line is given by a vector 
field 

m:B 2 ^ S 2 

defined in a circular cross-section Vt = B 2 of a thin film and satisfying: 

V-m'^OinB 2 and m'(x) = x^ on dB 2 . (25) 

(For the Bloch line in a thin cylinder, the magnetization is assumed to be invariant in the thickness 




Figure 3: Bloch line. 

direction of the film and the word "line" refers to the vertical direction.) Since the magnetization 
turns in-plane at the boundary of the disk B 2 (so, deg(m', dil) = 1), a localized region is created, 
that is the core of the Bloch line of size e, where the magnetization becomes perpendicular to the 
horizontal plane (see Figure [3]). The reduced energy (|22|) for a configuration (|251) writes as: 

E red (m) = [ \Vm\ 2 dx+\ [ m\dx. 
Jb 2 £ Jb 2 

The Bloch line corresponds to the minimizer of this energy under the constraint (|25p . Remark 
that the reduced energy E red controls the Ginzburg-Landau energy, i.e., 

/ g £ (m')dx < E red {m) 
Jb 2 
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since |Vm'| 2 < |V(m' \m^)\ 2 and (1 — |w'| 2 ) 2 = m\ < m§. Due to the similarity with the Ginzburg- 
Landau type functional, the Bloch line corresponds to the Ginzburg-Landau vortex and the ener- 
getic cost of a Bloch line (per unit- length) is given by ([6]) : 

£ rerf (Bloch line) = 0(| log e|) 

with the exact prefactor 2ir (see e.g. [5]). 

Micromagnetic boundary vortex. Next we address micromagnetic boundary vortices. A 
micromagnetic boundary vortex corresponds to an in-plane transition of the magnetization along 
the boundary from v to — v , see Figure |H The corresponding minimization problem is given by 



E red (m) 



|Vm| 2 dx 



log VI 



{m! ■ v) 2 dH 1 



J on 

within the set of in-plane magnetizations m : Q — > S . The minimizer of this energy is an harmonic 
vector field with values in S* 1 driven by a pair of boundary vortices. These have been analyzed in 
[Ml [T6] . The transition is regularized on the length scale of the exchange part of the energy, i.e. 



the core of the boundary vortex has length of size 
by 

^red 



'1 



I log *H 



The cost of such a transition is given 



E re (Micromagnetic boundary vortex) = 0( 



log- 



V 



log VI 



with exact prefactor ir. (Note that the boundary vortices in Theorem [3] correspond in fact to "half 
Bloch lines where the vector field is tangent at the boundary, i.e., m' ■ v = on dCl; therefore, 
their structure is different from the one of micromagnetic boundary vortices, but with the same 
energetic cost.) 





Figure 4: A micromagnetic boundary vortex 

Claim: In the regime (f!2U]). then 

either E red (Micromagnetic boundary vortex) < _E red (Neel wall) 
or E red (Micromagnetic boundary vortex) < E red (Bloch line). 

Indeed, assume by contradiction that the above statement fails. Then one has 



7?|l0g77| 



< 



log- 



n 



iog ; 



and 



log - ^ 

e 



los 



I log VI 



(26) 



(27) 



In the regime (f2T)]) . one has e 2 « e « r), therefore (|2"6"]l turns into 

"TT — | ^ log log-, 
77] log T7| e 
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while (|27|i implies that 

log - < log — . 

£ 77 

Now it is easy to see the incompatibility between the last two inequalities as e, 77 — > 0. 

Our toy problem: The model we presented in the introduction consists in considering config- 
urations without lateral surface charges, i.e., (JT]) holds true. In this case, our energy functional 
E e 2„{rri) coincides with the reduced thin- film energy E red since h ac (m) induces the stray field en- 
ergy (|23p as in ([2]). However, ([1]) would be physical relevant for a global minimizer only if boundary 
vortices were more expensive than both the Neel walls and Bloch line contribution. As explained in 
the above Claim, this assumption is violated in the regime (|20|) . Therefore, our energy functional 
is not adapted for studying global minimizers in the regime (|20j) . but rather for metastable states 
that satisfy |T]). 

Recently, the regime £ red (Micromagnetic boundary vortex) -C E red (Neel wall) -C E red (Bloch line) 
was investigated in Ignat & Kniipfer [10] for thin films of circular cross-section. It is stated that 
the global minimal configuration for that geometry is given by a 360°— Neel wall that concentrates 
around a radius so that it becomes a vortex (the Landau state of a disk) at the mesoscopic level. 



4 Some preliminaries 

The result stated in Theorem [2] is an extension to the S 2 — valued magnetizations of the following 
compactness result for S 1 — valued magnetizations obtained by the authors in |llj : 

Theorem 4 (Ignat & Otto UJjj) Let B n be the unit ball in R n , n = 2,3. For every small r\ > 0, 
let m' n : B 2 S 1 and h v = (h' n , /i 3jl; ) : B 3 ->• R 3 be related by 

J (h' v - VC + fca,^) dxdz = J (\7-m' v dx, VCeC c °°(B 3 ). 



Suppose that 



1 f ,0 , , C 



\V-m'\ 2 dx + - \hJ 2 dxdz< ~ , , (28) 

' vJbs ?y|log?7l 

for some fixed constant C > 0. Then {jn^j^o "is relatively compact in L 1 ^ 2 ) and any accumulation 
point m' : B 2 — > M. 2 satisfies 

\m'\ = 1 a.e. in B 2 and V • m! = distributionally in B 2 . 



In the proof of Theorem |3l we will use the following result due to Jerrard [12] for the concen- 
tration of the GL energy (JT} around vortices (see also Sandier [IS] , Lin JTS] ) : 

Theorem 5 (Jerrard [121) Let a € [0,1) and d > be a positive integer. There exists e = 
£o{d, a) > such that for every < e < eq, if m! : Vt — > R 2 satisfies the following conditions: 

\m'\ > i on {a; £ H : dist(a;,gO) < r*(e)} /or some r* (e) £ ( | i g £ |4 ' ^ 29 ^ 

|deg(m',an)| = d 
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and 

g e (m') dx < 2n(d + a) \ log e | , 



then there exist n points x\, . . . ,x n € Q with dist (xj , dQ) > r*(e), j = 1, ...,n and positive 
integers d\, . . . , d n > suc/i f/ia£ i/ie n 6aHs {_B(a;j,r*(e))}i<j< n are disjoint, 

n 

Y,dj=d 

3=1 

r*(e) 

g £ (m')dx > 2^1 log 1 - C{d,a), j = 1, . . . ,n, 

where C{d,a) is a constant only depending on d and a. 

In the above theorem, Q is any open bounded set (without any regularity condition imposed 
for the boundary dfl). This is due to hypothesis (|29p of having a security region around dfl. By 
degree of a C 1 — function v : C — > S 1 defined on a closed curve C with the unit tangential vector r, 
we mean the winding number 



deg(v,C) = — [ det(«, d^dH 1 . 
2tt Jc 



If m! : C ->• M 2 is a C^-function with \m'\ > on C, we set deg(m',C) := deg(-f^- 7T , C). The notion 

|m | 

of degree can be extended to continuous vector fields and more generally, VMO vector fields, in 
particular iJ 1 / 2 ^, 5* 1 ) maps (see Brezis & Nirenberg [2]). 



5 Proof of Theorem [3] and Proposition [T] 

First of all, let us define the security region around dfl together with some notations that we use 
in the sequel: 

Definition 1 Let fl is a simply- connected bounded domain of C ' boundary. The security region 
around <9f2 is the maximal set of points around dfl ( in the interior and outside £1 ) covered by 
the normal lines at dQ before any crossing occurs. We call depth of the security region to be the 
smallest distance to the boundary dQ where a crossing of two normal lines occurs and it will be 
denoted by R(dfl). 

Let R = R(dfl) be the depth of the security region around For r £ (0,i?), we denote the 
interior subdomain tt r C f2 at a distance r from the boundary, i.e., 

S! r = {i£fi : dist (x, dtt) > r} and 8Q r = {x e Q : dist (x, OQ) = r} (30) 

be the boundary of this subdomain. For r £ (— R, 0), we write dfl r to be the symmetry of 9f2_ r 
across the boundary dQ — dflo and £l r D Q be the extended domain surrounded by dfl r . 

Let I = - H 1 (9il) be the length of <9f2. Set w : [0, 1] — > dfl be a C parametrization of 9f2 such 
that |w(s)| = 1 with w(s) = ^(s) and let v(s) — w(s) 1 - be the outer unit normal vector on 9il 
at w(s). Since w(s) ~ ^r(s) is parallel to v(s) for a.e. s S [0, 1], we will always write 

w(s) = w(s)u(s) 
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where w(s) is the signed length of the vector w(s) with respect to u(s). Notice that |u)(s)| < g7g^ ■ 
In the security region around <9£7, a point x writes in the new coordinates as: 

x = F(s,t)=w(s)+tv(s), s G [Q,Z],t G (-R(dQ),R(dQ)). (31) 

Note that for interior points x G fi, the corresponding normal coordinate £ is negative. We define 
the symmetry transform $ in the security region around dtt: 

<f>(F(s,t)) = F(s,-t) se[0,l],te(-R(dQ),R(dn)). (32) 



A first ingredient that we need in the proof of Theorem |3] is a mirror- reflection extension across 
the boundary dfl. 



in '(<I>(x)) 




m'(x) 



Figure 5: Mirror- reflection extension. 



Lemma 1 Let Roo > 0. There exists eo = eo(R co ) > such that for every < e < Eq, the 
following holds: 

Let fl be a simply- connected bounded domain of C ' boundary with the depth of the security 
region R(dfl) > R^. Let $ be the symmetry transform across the boundary dQ defined in (|32p. 
In the security region, we consider the interior curve 

7 = an i 



(see notation (|30j) ) and m! : Q — > B 2 is a C 1 vector field that satisfies ([T]) 7 

|m'| > 1/2 on 7 and deg(m ,7) = 0. 



Then there exists an extension vector field ml : 0__ 



of ml into the extended domain 



D of boundary 



7 = $(7) = dtt 1 



such that 



rh = m! in f2, \rh'\ > 1/2 on 7 and deg(?ri',7) = 2 

1 



g e [m')m\y)- / g e {m') dU\x) 



< C \loge\H\dn) 



J 

7<I> 



<S>(W) 



g e (m')dy- / g e (m')dx 



IT 



jloge 
1 



g.im^dU^x)) (33) 



|loge| 



g s (m') dx 



IT 



(34) 
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where W C O \ 17 1 is any open subset of 17 and C — C(Roo) is a positive constant depending 

log i 

onZi/ on i?oo . 

Proof of Lemma Q] . We use the notations introduced at the beginning of this section. We have 
that \w(s)\ < R (g n j < Moreover, differentiating (j^TT) . we have that for a.e. s £ [0,/] and 

t€(-R(dtt),R(dn)), 

DF(s,t) = ( a s (t)w(s) ) and DF~ 1 (s, t) = ( ^jw(s) ) T , (35) 

where 

a s (t) := 1 — tw(s). 

By (j3"2")l and (|35l) . we compute that: 

SJt) := D$(x) = — ^— w(s) ® w(s) - Id for a.e. s € [0, 1] and t <= (-R(dn), R(dil)). (36) 

The matrix S s (t) is symmetric and its inverse is given by S^i) -1 = S s (—t). The mirror-reflection 
extension m! of m' is defined as: 

m'($(i)) := 5,(0)m'(i) = 2m'(i) • w(s)w(s) - m'(i) for x e 17 \ 17 fl(an) . (37) 

(We use that a<E>bc = (6 • c)a, for any a, 6, c € K 2 . ) Remark that the condition (JlJ implies that the 
mirror-reflection extension does not induce jumps at the boundary. Moreover, |m'($(x))| = |m'(x)| 
since S s (0) — 2w(s) <£> w(s) - Id is a reflection matrix (i.e., it is symmetric and orthogonal). 
Therefore, |m'| > 1/2 on 7. 

The goal is to estimate the energies / g E (rh')dy and / g e {rn) dT-L l . We start by computing 
the Dirichlet energy of the extension m! . For that, we differentiate (|3"T)) in the coordinates (s,i): 

D\m'($(z))j = S s (0)Dm'(x)DF(s,t)+2^ V(s)m'(x) ) , 

where 

V(s) := w(s) <g> w(s) + w(s) <8> w(s). (38) 
Since Z>( m'($(x)) ) = Dm'($(x))D$(x)DF (s, t), multiplying by DF(s, t)~ 1 S s (—t), it implies 



that 

M'($W) ^i 133 S s (0)Dm'(x)S s (-t) + —L-V(s)m'(x) ® w(s). 

a s (-t) 

Since 

(flm'($(i))) T = S* s (-i)L>m'(x) T 5 ;5 (0) + — — -w(s) ® V(s)m'(x), 

a s {-t) 

it follows that 

|Dm ($(x))| 2 = tr(£>m' ($(i))M' (*(x)) T ) 



tr(S s (0)Dm'(x)S s (-t) 2 Dm'(x) T S s (0)) + ^\V (s)m' {x)\ 2 



ot s (-t) 

4 -tr(S s (0)Dm'{x)S s {-t)w{s) ® 7(s)m'(x)) 



a«(-i) 

= 1 + 11 + III. (39) 
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For the first term in (l39l) . we compute that 



„ , , N? HQ 4:tw(s) ... ... 

S s {-t) 2 = ^w(s) ® w(s) + Jd. 

a a (-i) 

Since t^SAS -1 ) = tr(A) and tr(Av (g) Aw) = |Aw| 2 < |A| 2 |w| 2 for any two matrices A and 5 in 
R 2x2 with S invertible and any vector w € K 2 , we deduce that 

I<(l + ^m^)\Dm\x)\\ (40) 



a s {-ty 



For the second term in (j31?)) . we have that ^(s)! 2 = 2|w(s) ® w(s)| 2 = 2|u>(s)| 2 and therefore, 

II < -*\V{,)\*\w!{x)\* < ^*K(z)| 2 . (41) 
a s (-t) a a [—t) 



For the third term in (|39|) . we compute that 

&(-*)*(•) ^ 5§yw( S ) and S s (0)V(s)^Ms)(^ j). 

Using that ir(A& ® c) = c ■ Afo and tr(A) = tr(S s (0)AS s (0)) for any matrix A in R 2x2 and any 
vectors b, c € M 2 , we deduce that 

jjj = 4a '® frf-Dm'CaOwfo) ® (S,(0)7(«X(i)) 
a s (-i) V 

a s (-t) 

= _ ^(t)w(s) ml{x)± {Dml{x)Ms)) 
a s (-t) 
ia s (t)\w(s)\ 

< |m (a)] \Dm (x)\. (42) 

a s (-<) 



Since | det£>$(a;)| I 1 , we deduce by Jgg), gD}, gj) and @2), 



sit) 

|^(a>(,))| 2 | det(^ ( *))| < ^ f i + »M) |^ ( ,)| 2 + ^f^-K (,)l 2 

+ ^^|to'(x)|| j Dto'(x)|. (43) 
Therefore, for every open set W C f2 \ fl i , we obtain by Young's inequality, 

log £ I 

\Drh'(y)\ 2 dy = / | detD$(x)\ \Dm' (<S>(x))\ 2 dx 
<s>(w) Jw 

^ / i( 1+ n iOI^'O'OI + C|loge||m'(a;)| 2 \ dx 

Jw 



logs] 



< (1 + -nr^— p) / |-Dm'(g)| 2 da; + CW 1 (9ft), 
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with C = C(iZoo) > and e < e(i?oo)- (We use that H 2 (VF) < p^W 1 ^).) Also, 
/ (1- |m'(y)| 2 ) 2 dy = / | det D$(x)| (1 - K(z)| 2 ) 2 rfx 

J&(W) iff 

^ 1+ M^) / (1-K(a:)| 2 ) a dr. 
|loge| J w 

Therefore, we obtain: 



g E {m')dy<{l + -f—) f g B (m')dx + CH 1 (dil). 
I log e| ,/w 

By changing f to —t in the above argument, the inverse inequality also holds: 
g e ( m ')dx<(l + -^-) f g E {m')dy + Cn l {dn). 

I log £ I y*(w) 



(44) 



Thus, inequality (|34|l immediately follows. For proving inequality (|33|l . we proceed in the same 
way: Since — n g e i ) = n> — i lo g E i ^ is a Lipschitz parametrization of 7, we compute 



i-($(F( S) --^))) 
as I log £■ I 



l36l 



det D$(F(s, 



|loge 



))l 



d 1 



and we have by l|4"3] 



g s (m'(y))dH 1 (y) = / | det £>$(a;)| ffe (m'($(a:))) cW 1 (x) 
$( 7 ) J 7 

^(I + tt^-t) / 5£ (m'( a; ))dW 1 (x) + C|log£|^ 1 (9Q). 
I loge| J 1 



By symmetry, (1331) follows immediately. 

It remains to prove that if deg(m',7) = 0, then deg(m',7) = 2. For that let <po : [0,1] — > R be 
the lifting of w, i.e., w(s) = e llp °( s \ Obviously, 



deg(w) = ^(^ o (0-^(0)) = l. 

Z7T 



(45) 



On the curve 7, we know that m! £ C 1 (7,R 2 ) and we write m' = pv with p = \m!\ > 1/2 and 
v : 7 — > S 1 . Then p, v € C 1 (7). Then deg(i;, 7) = deg(m', 7) = 0. In this case, the theory of lifting- 
yields the existence of a lifting <p g C 1 (7,R) such that v = e lip . If f := i lo g e | > then F(-,— t) is a 
parametrization of 7 and we have 



= deg(«, 7 ) = —(tp(F(l, -t)) - ^(F(0, -*))). 
Notice that the reflection matrix S s (0) has the following form: 



(46) 



s /qx El j cos2<,5o(s) sin2<^o(s) j 
I sin2<po(s) — cos2(/?o(s) / 

That implies the following writing of m! on the curve 7 = $(7) parametrized by F(-,t): 



m'(F(s, t)) = p(F(s, -t))S.(0)v(F(s, -t)) - p(F(s, -t))e 



i 2 Vo (s)- V (F( S ,-t) 
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Therefore, by (|45|) and (|46|). we conclude that 



deg(m',7) = 2. 

□ 

We now prove the concentration of the Ginzburg-Landau energy on a small region (either one 
interior vortex, or two boundary vortices) under the condition ([!]) in the regime (1151) : 

Proof of Theorem \3[ Let R = R(dft) be the depth of the security region around dfl. We proceed 
in several steps: 

Step 1. Find a good set of boundaries. We define the set I of distances r € (e, R) such that we 
control the energy of m! on the boundary d£l r (and consequently, the modulus \m'\ via (09])), i.e., 



/= |r€ (0,R-e) : J g,{m l )dU 1 < |loge| 3 | 



(47) 



How large is the set 7? We show that for each interval Jc (0, R — e) of length £ ^> | l0 g E ^ , there 
exist infinitely many distances r belonging to J PI J. More precisely, we have for small e > that 



H\lnJ)>-. (48) 



Indeed, one has 

CD 



gu r r 

4tt| loge| > / / g e {m')dn 1 dr>\\oge\ 3 n 1 
Jj\i Jan r 



(J\I) 



which yields 'H 1 (J \ I) < | < | for small e > and therefore, (l48t holds. Moreover, remark 

that > i for every r € I, if e > is small enough. Indeed, since r < R — e it means that 
'H 1 {d^l r ) > rL 1 {dB{Q,e)) > e. Denoting by p := \m'\ and min := min{p(x) : x e 9f2 r }, it is easy 
to check that (see Lemma 2.3. in [12)): 

|loge| 3 > / g^m^d-H 1 > f ( \d T p\ 2 + -1(1 - p 2 ) 2 ] dU 1 > -(1 - mm) 2 , (49) 

where r is the tangent unit vector at dil r . Thus, one concludes that (1 — min) 2 < e\ loge| 3 -C 1, 
i.e., min > 1/2 for small e > 0. (Relation (j49|) is obvious if p is constant (equal with min). 
Otherwise, the GL energy of the modulus p controls the following quantity ^ J/ m (p)(l — V 2 ) dy on 
the image set Im(p) of p and forces p to take values close to 1.) 

Bv (|48]) . we can choose ri € I fl (^n^r j i i g E i )• W.l.o.g., we may suppose that 

1 

Tl I logej " 

We distinguish two cases in function of deg(m', dfl ri ). 

Step 2. FFe assume that | deg(m', 9f2 ri )| > 0. W.l.o.g. we may suppose that d := deg(m', dSl ri ) > 
1. 
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Extension. We will extend the vector field ml 
fl ri _ r * D f2 ri with 



by a vector field m! defined on the larger domain 



such that m! — m! in Q ri and 



deg(m', d^n-r') =d>l, 
\m'\ > ^ in O ri _ r » \ Q ri 

[ g e (rh')dx < 2tt(1 + a)|loge| + C. (50) 

For that, using the notation (|3Tj). for each point a; = F(s,t) £ fl ri _ r * \ S! ri (here, i < 0), we 
consider y x = F(s, — ri) G <9i7 ri to be the normal projection of x on dft ri with dist (a;, dil ri ) = 
\x — y x \ = \t + ri|. Then we define m' : fi ri _ r * — > R 2 as m! — m! in 51 ri and 

m! \x) := m!{y x ), for every x £ Sl ri _ r » \ Q ri . (51) 

Since ri £ /, (|49|) implies |m'| > 1/2 on f2 ri _ r * \ Q, ri , deg(m' , dfl ri - r *) = d > 1 and 

/ g £ (m')dx = [ [ g e {m!)dU 1 dr 

< CV* / g^m^dU 1 

Jdfl ri 

< Cr* | log£| 3 = C. (52) 

Thus, by (|T5|) . we obtain (|50l) . By Theorem [5] and (|50[) . we deduce that d = 1 and there exists a 
point x\ £ f2 ri such that 



/ g e (m')dx>2n\log-\-C(a), 



where C(a) is a constant depending only on a. Therefore, we obtain via (1521) that 

/ g e {m')dx >2%\\og — \-C(a). 

JB(x 1 ,r*) £ 

Since B(xi,r*) C f2, the conclusion (|16[) follows. 

Step 3. We now dea? wii/i i/ie other case deg(m', <9Sl ri ) = 0. 

Mirror-reflection extension. We consider the symmetry transform $ defined in Lemma [T] across 
the boundary dtt together with the mirror- reflection extension m! : Q r2 — > B 2 where T2 = — r\. 
Then Lemma [T] yields 

/ g e {m') dx < 2tt(2 + 2a) | log e\ +C(dQ), (53) 
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\rh'\ > i on 90 r2 , the degree of m' on the boundary <90 r2 is equal to 2 and 



/ g e {m')dn 1 < 2\\oge\ 



for e small enough. The extension argument in Step 2 leads via Theorem [5] to the concentration of 
the Ginzburg-Landau energy of ml into two vortex balls 8(2:2, r*) and 8(2:3, r*) with 2:2, 23 S O r2 
and there exist two non-negative numbers > d% > 0, di + dj, = 2 such that 

/ g e (m')dx>2ird j \]og — \-C, J=2,3. (54) 

JB( X] ,r») £ 

(The assumption a < 1/2 is needed so that 2 + 2a < 3.) 

Case 1: di = 2 (i.e., there is one vortex ball of degree 2 in 0,. 2 ). The level of energy (fT5"|) rules out 
that 8(x 2 , r*) C O. By Lemma[IJ it also means that B(x2, r*) is not included in O r2 \0 (otherwise, 
the symmetry of the energy distribution around the boundary would imply again that the reflected 
domain Q(B(x2, r*)) charges the energy more than the level fTS")) in the interior of O). Therefore, 
B(a;2,r*)n90 ^ 0. Choose 24 = a;| € B{x 2 , r*) n 90. Then B(x 2 , r*), $(B(x 2 , r*)) C £(24, 10r*) 
and by Lemma Q] and (l54|) , we conclude that 

g s (m') dx > — [ / g e (m / )dx + I g e (fh')dx\ 

B(x^,Wr*)nn 2 yjs(2;2,»-*)nn J<s>(B(x 2 ,r*)\n) J 

g E (m')dx + / 5 £ (m')c?x — C 

B(x 2 ,r*)nn J B(x 2 ,r*)\n J 




g e (m')dx-C > 2ir\ log — | - C*. 

lB(x 2 ,r*) e 

Here, <P77|> holds and -B^l, 10r*) = 8(2:3, 10r*) are the boundary vortex balls. 

Case 2: d2 = d$ = 1 (i.e., there are two disjoint vortex balls of degree 1 in r2 ). If (Til))) 
holds, then we are done. Suppose that (p~6|) is not satisfied. Then we want to prove (| 1 T|) . As 
in Case 1, the symmetry of the energy distribution around the boundary implies via Lemma U 
that none of the balls B(x2,r*) and -8(2:3, r*) is included in O or O r2 \ O (otherwise, p^]) would 
hold). Therefore, B{ Xj ,r*) n 90 ^ for j = 2, 3. Choose x* e B{ Xj ,r*) n 90 for j = 2, 3. Then 
B(a;j,r*),$(B(afj-,r*)) C 8(2;*, 10r*) for j = 2,3. As before, by LemmaQ]and <[53]l. we conclude 
that 

/ g e {m')dx>\( I 9e(m')+ 

J (B(x*,10r*)uB(xl,10r*))nn L \ J [B(x 2 ,r* )UB(x 3 ,r*) J nn 



(Here, we used that $ ^8(2:2, r*) \ Oj n <f> I -8(2:3, r*) \ OJ =0 since the two balls B(x 2 , r*) and 
.8(2:3, r*) are disjoint and lie in the security region of 90.) □ 



*(B(:E2,r*)\n)u4(.B(x3,r*)\n) 

>\ I g £ (rh')dx-C 

Z Jfl(i2,r*)UB(i3,r*) 

> 2tt| log — I -C. 
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The natural question is whether the lower bound for the energy of a boundary vortex given 
in Theorem [3] is optimal. A positive answer is supported by the following result: we prove that 
the loss of energy of order i g°f g 6 | (with respect to n\ loge| which is the exact half energy of an 
interior vortex) may be achieved for C 1 domains. 

Proof of Proposition^} The aim is to construct a boundary vortex on <9f2 centered at the origin. 
Step 1. Construction of m' e . We define the two-dimensional vector field m! e that is tangent at dQ 
and its phase ip £ is linear on every arc of circle {\x\ — r} n £1 with r G (0, 1/200). More precisely, 
let 

f e i Ve {x) if £ e ft and e < |x| < 1/200, 

if x £ il and < Id < e, 



m' £ (x) 



where the phase tp e is given in the polar coordinates as follows: tp e (r,-) : (—j(f),j(r)) 
(— 7r/2,7r/2) is an odd function (i.e., tp e (r,Q) = ~ip e (r,—9)) and it is linear in 9, 



(r) 



7(r) 



The phase correction 86 : (0, 1/200) 



JS6(r 



\ 1 + irj'(r) 



\/l + (r 7 '(r)) 2 



for every 6 G (-j(r), j(r)), r G (0,1/200). 
( — f-,0) due to the condition (fTJ) is defined as 
for every re (0,1/200). 



(Indeed, one can easily check that m' e is tangent at dft.) 

Step 2. Estimate of the GL energy outside the core region. We first estimate the energy of m' e 
away from the core, i.e., V\ = {x e ft : e < \x\ < 1/200}. For that, we need the following 
computations: for r 6 (0, 1/200), 



7 '(r) 



■log (log log i) 2 ' 



89(r) = — arccos ■ 



t55l 

\89(r)\ < 2\rj'(r)\ < 



log ~ (log log i) 2 ' 



(55) 



(56) 



y/i + (n'(r)y 

where we used that arccos : [—1, 1] — > [0,7r] satisfies yjl — t 2 < arccosi < 2^1 — t 2 for t G [i, 1]. 
By a change of variable r — e s and differentiating (|5fc>[) in the new logarithmic variable s, we deduce 
that 



ar 



as 



1 



f7(e s ) 



d 2 , s, 



(57) 
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and we have for a universal constant C > 0, 
/ 9s(m' £ )dx = / \Vif E \ 2 dx 

JVi JVi 



1/200 n (r) /la |2 
-7(r) 



< 



r 1/20 ° f d } ch 

< J |2 7 (r) + C 7 (r)r 2 ( — ^(r)) 2 j — + C 



r 



s=log r, 113 /"log 5SJ f C / d 2 , , 

< / |2 7 (e s )+7(e s )- (_ 7 ( e «) ) ^ + C 

/ (•< -r—H 1 nrl— m< if K ( 



logs 



log \s\ J \ \ds 2 



<n\\oge\-2 —- 1- -)ds + C 

log | log £ | 

(Here we used that f* Q ds ~ as i -> oo.) 

Step 3. Estimate of the GL energy inside the core region. Now we estimate the energy of m' E on 
the core, i.e., T>% — {x 6 £1 : < \x\ < e}. Using the same argument as above and the change of 
coordinates r = e s , we compute 



Vm'l 2 dx= I I | ' ' T' + ) rdddr 



T> 2 JO J-y(r) 



r 2 |V<^| 2 , 1 



E 2 



dHJ,4sSJ f e / Cr r 3 Tfr) / d \ 2 \ (13 CV 



and 



/ l(i-K| 2 ) 2 ^ = o(i). 



Step 4. Conclusion. The -ff 1 vector field satisfies the required properties in Proposition [T] 
However, m^. is not C 1 . In order to construct a smooth m^, we define / : M. [0, 1] a smooth 
function such that f(t) — if t < and f(t) — 1 if t > 1. Now, it is enough to change the vector 
field defined at Step 1 only in the core region as follows: m' £ (x) = /(HyM*) if Q < |a;| < e . 
□ 



Remark 2 The GL energy of a boundary vortex placed in a corner is proportional with the 
corner angle. Therefore, the loss of energy of order lo g°f tg| f° r C 1 boundaries (see Proposition 
[1} increases to an order of |loge| for Lipschitz boundaries. More precisely, let f2 = {(xi,x%) : 
x\ € (—1, 1), x 2 > \x\ \ tan where a G (0, n) is the corner angle of f2 at the origin. For every 



22 



< e < 1, we consider the following approximation of a vortex: 

m' e {x) -- 



jfy ife<|a;|<l, 



if < |as| < £. 
Then m' £ satisfies (P) on dfl n B 2 and 

/ g e (m' e )dx < a|loge| + 0(1). 
JnnB 2 

6 Proof of Theorem [2] 

We will work at the level of sequences of parameters Ek and rjk and a sequence of magnetizations 
rrik satisfying the assumptions in Theorem[2l We will prove the Theorem in a slightly larger regime 
than (fT()|) ; more precisely, it is enough to assume that 



4 £ V, (58) 

for some constant (3 G (0,1) such that 

P < 1 - a. 

By H3D, let A > be such that 

.4 

Vk | log T] k 

By (fTTj) , there exists a > such that 



£e*,7» ( m fc) - 2tt| loge fe | < - for every fcelf. (59) 



a log | log e fc | < — (60) 

W log i?k I 

for every k G N. Let [/& : R 3 — >• K be the stray field potential associated to mk defined by (|24|) for 



(V, -§-)Uk that satisfies 



/ f|VC/ fc | 2 + |^| 2 ") dxdz^i / iVr^^.m'J 2 ^. (61) 

JR 2 xR V OZ J Z J R2 

By the Lax-Milgram theorem, the potential Uk exists and is unique in the Beppo-Levi space (see 
Dautray and Lions [3]): 

BC = {U : R 3 -> R : (V, %-)U G L 2 (R 3 ), — ^— r G L 2 (R 3 )}. 
L v <9z ; v ; 1 + |a:| 

We proceed in several steps: 
Step 1. Location of the vortex balls of m! k . Let r* = 1 / 1 log £fe | 3 . By Theorem[3l there exist at 
most two points Xk ,4 6 such that 

f g Ek (m' k ) dx > 2n\\og— | - C7 > 2tt| loge fc | - lOOlog | loge fc |, (62) 

for k sufficiently large. Obviously, up to a subsequence, {xk},{ik} C fl converge to two points 
Xo-, Xo G f2 and we have that for every small a > 0, 

B(x k ,r*) C S(x ,(t), B(x k ,r*) C S(x ,cr) 
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for k sufficiently large. 

The set V = B{xq,<j) U B(xo,a) is the location of the essential topological defects of each 
m' k . Now the goal is to prove that {m k } is relatively compact in \ T>). The idea is to 

approximate m' k away from T> by S 1 — valued vector fields, denoted by ML that satisfy the hypothesis 
of Theorem 2J For that, let B C ft \ T> be an arbitrary ball. To simplify the notation, let 
B = B(0, 2) C M 2 be the ball of radius 2 centered in the origin. Since m 2 k > k = (1 — Im^ 2 ) 2 , 
the energy level on B is bounded as follows: 

/ \\7m k \ 2 dx+\ ( (l~\m' k \ 2 ) 2 dx+- [ \\V\-^ 2 (V ■ m' k )\ 2 dx 
Jb £ k Jb Vk JR 2 

< E e k ,r, k (m k ) - / g £k {m' k ) dx 
Jvnn 

{T2j,{59},!62} a 

< min{ — -,27ra|loge fc |} + 100 log | loge fc | 

Vk\ log Vk\ 

< min^— f ,27ra|log£ fe | + 1001og|loge fe | \ (63) 

[r] k \logr] k \ J 

for k sufficiently large and A = A(a + 100)/a (by ([6TI)l V 

Step 2. Construction of a square grid. For each shift t £ [0,e^), write 




Figure 6: The net of horizontal lines. 

V t := {(x 1 ,x 2 ) £ B : x 2 = t ( mod e k )} 

for the net of horizontal lines at a distance e k in B. By the mean value theorem, there exists 
t k G (0,£ k ) such that 

/ 9e k (m k ) dU 1 < \\ g £k (m' k )dx. 
Jv tk e k Jb 

If one repeats the above argument for the net of vertical lines at a distance e k in B, we get a square 
grid lZ k of size e k such that the convex hull of lZ k covers the unit ball B 2 C B and 

f feK)^' 1 ? -H-ri?— 7.3^)- m 

By the same argument as in (|49|) . the estimate (|64|) together with /3 < 1 implies that lZ k C {Ito^I > 
1/2} for k large enough. 

Step 3. Vanishing degree on the cells of the grid. In order to approximate m' k in B 2 by S 1 — valued 
vector fields with uniformly bounded H 1 — norm, it is necessary for m' k to have zero degree on each 
cell of the square grid lZ k . This property of vanishing degree is shown in the following lemma: 
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Lemma 2 Let < a < 1, < /3 < 1 — a and C > 0. There exists Eo = Eo(a, (3, C) > such that 

2 ' 2 



for every e G (0,£q) the following holds: if Z = (— V) 2 * s the cell of length and m! : Z — > £? 2 



is a C 1 vector field such that 

[ gJm^dH 1 < C \ lo % £ \ and [ gJm') dx < 2tto\ logd + (71og I logel, 
i/ien deg(m', 9Z) = 0. 

Proof of Lemma [2] . Note that the same argument as in (H^l) implies that |m'| > 1/2 on dZ, so 
that it makes sense to speak about the degree of m' on dZ. Note also that the quantity Clog | loge| 
in the upper bound of the GL energy on Z can be absorbed by the leading order term 2ira\ log e\ for 
a slightly bigger a > a so that the inequality f3 < 1 — a still holds. Therefore, we omit that second 
leading order term in the following. The idea of the proof consists in a rescaling and extension 
argument so that the imposed upper bounds on the GL energy rule out the existence of a vortex 
in the interior. Indeed, assume by contradiction that | deg(m', dZ)\ > (i.e., a vortex exists in the 
interior). By a change of scale, we define m! on the rescaled cell Zi/ 2 = (—1/2, 1/2) 2 as: 

m'(x) = m'(e x) if x G Z 1/2 = (-1/2, 1/2) 2 

and then, we extend m' to the larger cell Z\ = (—A, A) 2 with A > 1/2 (to be chosen later) as 
follows: 

m' (x) = rh' (y) if x G Z x \Zi/ 2 , 

where y € dZ 1 / 2 with y = tx for some t € (0, 1) (i.e., y is the closest point to x on the boundary 
dZi/2 that has the same direction as x). Therefore, |m'| > 1/2 in Z\\Zi/ 2 and | deg(m', dZ\)\ > 0. 
Letting S — e 1_/3 , we have that 



</,,( m') <lr j g e (m') dx < l °S 6 \ 



and 



g s {fh!) dx < C(X - h [ g 5 {fh') dH 1 

= CeP(X-h f geirn^dH 1 
1 Jdz 

<CC(X- 1 -)\loge\= CC ^_y ) \logS\, 

for some universal constant C > and a small e > 0. We choose A > 1/2 such that 

K:=^ + dCiX -^ <l 
1-/3 2tt(1-/3) 

(this is possible since by hypothesis, f3 < 1 — a). By summing over the above energy estimates, we 
obtain that 

/ g s {m')dx <2nK\\og5\. (65) 
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Since K < 1 , Theorem [5] implies the existence of a ball B C Z\ of radius A — 1/2 with 

gs(m')dx > 2tt| log<5| - C 



for <5 sufficiently small, which is a contradiction with (|65|) . □ 

As a consequence of Lemma [21 we deduce by (|63|) and (|64|) that our choice /3 < 1 — a implies 
that m k has vanishing degree on every cell of the grid lZ k . 
Step 4. Construction of an approximating sequence. We denote 

Pk = \m' k \ and m' k = p k v k . 

By Step 3, we can smoothly lift m' k on the grid, i.e., 

v k = ^ = e^ fc on ft fc and ip k € C 1 , R) . 

On each cell Z fe of length e k of the grid, we define 

M' h = e l * fc in Z k 

where $> k is the harmonic extension of <p k inside Z k , i.e., 

f A$ fc = in Z k 
\*k = <Pk ondZ k . 

Since ^/c can be smoothly extended around dZ k (because m' k / p k has a C 1 lifting around dZ k ), we 
deduce that & k g C 1 (Z k ). Note that the following inequality holds: 

\V<S> k \ 2 dx<Ce' 3 k ( \V Vk \ 2 dU\ (66) 

Z k JdZ k 

Indeed, after rescaling by e@ , we show the inequality in the unit cell Z\ — (— 1, l) 2 for the harmonic 
function $ in Z\ with the trace ip on dZ\ . We can assume that J dz f dH 1 — (otherwise, consider 
ip — j- dz ^ipd'H 1 ). For that, we consider a smooth cut-off function "J : [0, 1] — > R such that vE'(i) = 
for t < 1/2 and *(1) = 1 and the following extension § ext of ip in Zf. <& ext {t ■ x) = &(t)tp(x) for 
every t € (0, 1) and x € dZ\. By Poincare's inequality, one has 

f \V$\ 2 dx< f \V<S> ext \ 2 dx<C f {\Vip\ 2 + (p 2 )dH 1 <C f iVtpfdH 1 . 

JZ X JZi JdZ x JdZx 

The goal is to prove that the sequence {M' k } approximates {m' k } in L 2 (B 2 , R 2 ) and M' k satisfies 
([28]) for some associated stray field h k defined in B 3 . 

Step 5. Estimate \\V(M' k — mJJH^. Denoting by C a generic universal constant, we have 

\\7M k \ 2 dx= I \\7<$> k \ 2 dx 

z k Jz k 

< c4 / |Vp fcr W 

JdZ k 

= C4 f iVvrfdH 1 

JdZ k 

<Cef/ p 2 k \\7v k \ 2 dH 1 < Ce? [ \Vm' k \ 2 dH 1 (67) 
./az fc Jaz fc 
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since p k > 1/2 on lZ k . Summing up after all cells Z k of 1Z k , since the convex hull of IZk covers 
B 2 , we obtain by ([Ml) . 



\VM' k \ 2 dx<Cel 



g eh {m' k )<M}< 



C 



Vk\ log 7]k\ 



Combining with (|63[) . it yields 



\V(M' k -m' k )\ 2 dx< ° 

B 2 77fe | log 77fc I 



Step 6. Estimate \\M' k — mjjj^. By Poincare's inequality, we have for each cell Z k of 7Z k : 



z k 



ML - + M' h 

dZ k 



dx < Cef 



\WM' h \ 2 dx < Cef 



,.)Z k 



and 



z k 



m h - 



dZ k 



dx < Ce 
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z k 



\Vm' k \ 2 dx. 



(68) 



(69) 



(70) 



(71) 



Since v k = Ml on dZ k , by Jensen's inequality, we also compute 



ldZ k 









2 


m 'k) 


dx = / 


f {vk - m' k ) 


dx 




Jz k 


JdZ k 





< Ce 



2/3 



i)Z k 



(1-pkfdU 1 



J'dZ k 



<Ce 



13+2 



dZ k 



g^m'^dU 1 . 



Summing up (JTUJ) , CD} and over all the cells Z k of the grid TZ k , by 
that 



Cf 

\M' k ~m' k \ 2 dx< U£k 



2fl 



(72) 

and (|64p . we obtain 
(73) 

T]k\l0gT)k\' 

Step 7. Construction of an appropriate stray field h k associated to M' k in B 3 such that (|28[) holds 
for the couple (M' k , h k ). The choice of the stray field h k has the form 

h k := (V,^-)(U k + U k ) 
oz 

where U k is the stray field potential associated to m' k by (|6 1 [) and we consider U k E Hq(B 3 ) to be 
the unique solution of the variational problem 



B 3 \ Oz oz 



dxdz = I CV • (M k - m' k ) dx, VC € H^(B 3 ^ 
Jb 2 



(74) 



(It is a direct consequence of Lax-Milgram's Theorem in Hq(B 3 ).) Note that h k is indeed a stray 
field associated to M' k : B 2 — > S 1 on the unit ball. In order to estimate J B3 \h k \ 2 dxdz, we observe 
that 



d 



B :i 



dz' 



8 



\(y,—)U k \ 2 dxdz < \{V , —)U k \ 2 dxdz < 



dz 



C 

I log m\ 



(75) 



27 



by (|61[) and (|63p . It remains to estimate f„ 3 |(V , -§^)Uk\ 2 dxdz. For that, one should use an 
interpolation argument via and (I75|) . For that, we extend Uk by outside B 3 , so that the 

extended function (still denoted by Uk) belongs to _ff 1 (R 3 ) and the trace Uk € i/ 1 / 2 (R 2 ). 

R 2 

Moreover, we have 

2\ 



f \MV*U k \ 2 dx< \ [ l\vu k \ 2 

JR 2 1 JB 3 \ 



dU, 



dz 



dxdz = — 
2 



|VC/ fc | 



dU k 



dxdz. (76) 



Let us denote by T a linear continuous extension operator: 

T : H S (B 2 ) -> H S (R 2 ), a = 0,1. 
Then by interpolation, it follows that 



/ \M^T(M' k ^m' k )\ 2 dx< ( [ 
Jr 2 \ Js. 2 



\T(M' k -m' k )\ 2 dx 



1/2 



\VT{M' k -m' k )\ 2 dx 



1/2 



< c 



\ML-m' h \ 2 dx 



B 2 



1/2 



B 2 



\V{M' k -m' k )\ 2 dx 



1/2 



Combining with ([76|) . the choice C :— Uk in (f74|) yields 



B< 



dz 



dxdz — 



B 2 



< 



C/ fe V • (M' k - m' fc ) ote 
L> fc V • T(M' k - m' fe ) da; 

1/2 



< c 



j \\V\V*U k \*dx) (7 \\V\^T{M' k -m' k )\*dx 

JR 2 J \JR 2 

\vu k \ 2 



1/2 



dU k 



1/2 



B 2 



\M' k -m' k \ 2 dx 



dxdz 

1/4 



\V{M' k -m' k )\ 2 dx 



1/4 



Hence, 



d 



(55}, GU Cef El 1 C 



l(V,^)E4| 2 ^dz <" 
S3 ?7 fe | log rj k | 



< 



logr?fc| 



for k sufficiently large. Therefore, by ([75)1 and (|77l) . we conclude 

C 

|/i fc | 2 dzdz < -. 

B3 I log 77fc | 



(77) 



(78) 



By and (|78p. condition (|2"5)) is satisfied for Mi and the stray fields h k . Then Theorem 2] applies 
and implies that {M' k } is relatively compact in L X {B 2 ). Therefore, from (f73| . it follows that {m' k } 
also is relatively compact in L l {B 2 ). Since the ball B was arbitrary chosen in the complementary 
of V and we proved the relatively compactness result in the reduced ball B 2 = ^B, by a diagonal 
argument, we deduce that {m' k } converges in L l {Q. \ V) up to a subsequence. Letting now a — » 0, 
the conclusion of Theorem [2] follows. □ 
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7 Upper bound for the Landau state 



In this section we prove the upper bound stated in Theorem [T] for a stadium domain: 

Proof of Theorem [T] . The construction is carried out in several steps: 
Step 1: A Neel wall approximation. Let 

A := 77I \ogf]\. 

The parameter A corresponds to the core size of a 180° wall transition. More precisely, we consider 
the following Id transition layer (u\,V\) :K4S' that approximates a 180° Neel wall centered at 
the origin (see Figure [7]): 



ux{t) 



llo *f^ 1 if i*i < vr^y, 

I log A| I I — v j 

elsewhere, 



-Ul(t) if t < 0, 

yi-Hj(t) if*>o. 




OUi 



Figure 7: A 180° Neel wall approximation. 

The exchange energy corresponding to this transition layer estimates as follows (see DeSimone, 
Kniipfer and Otto [5] or Ignat [8]): 

1 dt 1 1 dt 1 " J R dt 1 VA|logA|y v ; 

In order to estimate the stray-field energy of the transition layer, let U\ be the radial extension of 
u x in R 2 : 



U\{xi,x 2 ) = u\{\Jx\ + x\ ). 
By if 1 / 2 (R)— trace estimate of an ij 1 (R 2 )— function, it follows (see details in [8] or (j96|) below): 

IMI W) ^ \ f. \™»\ 2 d X <,[\\ d ^\ 2 dr= I±£gl. (80) 



2 J R 2 J 1 dr 1 |logr?| 

We will construct a continuous vector field m : fi — >• S* 2 such that the upper bound in TheoremQ] 
holds and 

m'{x) = v^{x), m 3 (x) = if ieffi, (81) 
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where v is the outer unit normal vector on 90. Moreover, the function to will satisfy the following 
symmetry properties: 

m (x) = — m (— x), mi(x) = —mi(xi, —X2), 7712(21) = —mz^—xt, X2), ie!l. 



Step 2: Construction in f2j (the sub-domain defined in Theorem^) . We distinguish two regions 
in Hi (see Figured]): 

1 

Oi.i = {x e Cli : xi > 1 + 5} and £!i, 2 = {x e Oi : 1 < xi < 1 + 6} with <5 = 



log 77I 3 / 2 




Figure 8: The region fix. 
In we define to with values in S" 1 that behaves like a vortex centered in A = (1,0): 



m (x) = 



A 



\x-A\ 



, 7713(2;) =0 in 



By setting to' to be a 180° transition wall on dQ.2 H 50^2 (as in Step 1), i.e., 

m'(l,x 2 ) = {u x {x 2 ),v x {x 2 )) 1 - =e i9 ^\ to 3 =0 if x 2 & (-1,1), 

the vector field to is completely defined on dfli^ (together with (|8ip). Here, 0a is the angle 
transition between [0, 7r] of the 180° wall on dfli.2 fl <90 2 , i.e., #a(x 2 ) = and #a(— x 2 ) = tt if 
£2 € [— 1, —y/1 — A 2 ], and 



6\(x2) = arcsin 



|logA| l0g 



L=V e x {-x 2 ) = n - 9 x {x 2 ) if x 2 e [-\/l-A 2 ,0]. 
+ A 2 / 



Therefore, we define to' = e ltp , 7713 = inside O1.2 by a phase if that is uniquely determined by the 
boundary conditions on dil\ t 2 as an affine continuous function in x±: 



(1 + tyT- xl,x 2 ) =t arcsin ^1 - x\ + (1 - i)6 x (x 2 ), f € (0,1), x 2 G (-l,-\/l - <5 2 ), 



y>(l + (5/j, X2) = f arcsin ■ 



tp(x 1 ,X 2 ) = TT - tp(xi,-X 2 ), I € !ll,2, X 2 > 0. 



+ (l - t)e A (i 2 ), i G (0,1), x 2 G (-\/l-<5 2 ,0), 



We will denote by 



ots{x 2 ) = arcsin 
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the phase of the vortex at Oil 1,1 n dfli^- 

Step 3: Estimate of the exchange energy in Q±. First, we have: 



f |Vm| 2 da; = 0(|log5|) = 0(log|log?7|), (82) 
Jn 1A 



and 



|Vm| 2 dx = 2 J J \Vcp\ 2 dx = 0(6) 



\Vm\ 2 dx = 2 / / \W(p\ 2 dx 



|x 2 |<VT=^ x 2 e(-VT^,o) 



= 2/7° f tM^) - 9 x (x 2 )) 2 + 6(t^ + (1 - i)^) 2 ) dtdx 2 
Jo J-^T^W \<> dx 2 dx 2 



< 4 , , 

Introducing the notation a(^j-) = as{x 2 ), we compute: 

-VT^ 1 rl/S r l/8 



das 2 




' dx 2 ' 


1 02:2 / 



fVl-9 - ^ /-I/O /■I/O J 

/ W s {x 2 )dx 2 < a 2 {s)ds<A -tfa = 0(l) (83) 

Jo Jo Jo s + 1 



(where we use that arcsina; < 2x if x G (0, 1)), 







^(^rf^O ^, * 2 )=o(l), (84) 



<S AV zy ' v 5| log Al- 



and 



*l£l"" s j£ I^K^-i i^if^ ' 1 ' 

-Vi-^ ^/). „ ,j 1 



f V § id0x 
Jo dx 2 



2 dx 2 = 0(——) = o(- 



Therefore, 



A|logA| 77] log?7j 
/ |Vm| 2 ^ = o(— (85) 

Step 4: Construction in We define m by imposing the symmetry m(x) = —m{—x) for 2; € ^3. 
Therefore, by ^B2J) and ([55]). we have 

/ \Vm\ 2 dx= [ \\7m\ 2 dx = o{— i — r ). (86) 
■/n 3 Jn a 77] log r?j 

Step 5: Construction in £l 2 . We distinguish two regions in f2 2 (see Figure |9]): 

^2,1 = {x E f2 2 : |xi| e (2<5, 1)} and fi 2: 2 = {a € fii : < 26}. 
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0)3 



Figure 9: The region ^2,2- 

In 02,i, we define m with values in S 1 that behaves like a 180° Neel wall (as in Step I): 

m'(x) = (u x (x 2 ),v x (x 2 )) 1 - =e ie ^ x '\ m 3 (x) = for 2^ G (2,5,1), x 2 € (-1,1), 
m(x) = —m(—x) for x\ G (—1, —26), x 2 G (— 1, 1). 

Denoting by B r the disc centered at the origin of radius r, we decompose the domain 

^2,2 =B e UwiUW2UW3 

with 

lo 1 = {x&B 1 : |xi| < <5}\5 £ , 
cj 2 = {x G 2 \B 1 : \x!\< 5}, 
lo 3 = {5< \ Xl \ < 26} x (-1,1). 

In B e (the core of the vortex), we define 

m'(x) = sin( — |a;|) yj^J > fn^x) = \/l - \m'\ 2 (x) for x G B £ . 

In lji , we define m with values in S 11 that corresponds to the vortex away from the core: 

m!{x) = ' m 3( x ) — f° r x € B\ \ B e and |xi| < <5. 

In W2, we define m with values in S : m! — e tv ,m 3 — inside u> 2 . The phase <p is given as an 
affinc continuous function in x 2 determined by the values on the boundary du> 2 : 

ip(xi,— (1 — t) — t\J 1 — x\) = iarcsinxi, t G (0, l),£i G (0,(5), 

<p(xi,x 2 ) = 7r - <p(xi, -x 2 ), a; G w 2 ,xi G (0, <5),x 2 > 0, 

<p(x) = n + ip(—x), X G UJ 2 , X\ G (-5, 0). 

In W3, we also define m with values in S 1 where the phase (p is an affine continuous function in x\ 
determined by the boundary conditions on duj 2 : 

^(5 + 8t,x 2 ) = (1 ~ arcsin^ + t9 x {x 2 ), t G (0, 1), x 2 G (-1, - v 7 ! - <5 2 ), 

1 — VI — 

^(J + (5t,a;2) = {1 - t)a 5 (x 2 ) + t9 x {x 2 ), t€ (0,1), x 2 G (-^1 - * 2 ,0), 
^(xi,x 2 ) = 7r - -x 2 ), x\ G (<5, 2<5),x 2 G (0, 1), 

ip(x) = 7r + (p(-x), ari G (-2(5, -<5), x 2 G (-1, 1). 
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Step 6: Estimate of the exchange energy in £^2 - We start by estimating the exchange energy in 
Sl2,i and then, in ^2,2- By (|75|) . we have that 



Vm| 2 dx = 2(1 - 28) I ( |^| 2 + |^r )olt = o 



'n 2 ,i 

In ^2,2, we first have 
Then 



dt 



dt 



1 



7 ?l log 77I 



/ \\7m\ 2 dx = O(l). 
I |Vm| 2 da; = 27r|loge|-0(|logJ|), 



Wm\ 2 dx = 



|V^| 2 da; = 0(8) 



and 



|Vm| 2 dx = 



\V<p\ 2 dx^o( 



1 



7y|log?7 



)• 



(87) 

(88) 

(89) 
(90) 

(91) 



By ©, ([87 

of to: 



891) , (I90| and (1911) , we deduce the following estimate of the exchange energy 

1 



/ 



Vm| dx = 2tt\ logs| + o 



77I log 77I 



(92) 



Step 7: Symmetries of the stray field. Now we estimate the stray field energy of m. For that, let 
U € be the stray field potential in the Beppo-Levi space associated to m defined by (f24|) for 
(V, -§-)U that satisfies 

/ (\ WU \ 2 + \ir\ 2 ) dxdz= f U(x,0)V ■m'(x)dx = l [ |V|- 1/2 (V • m') % dx. (93) 
Je 2 xk V oz J Jn 2 J R 2 

Moreover, the stray field potential verifies: 



AU = 




if z ± 0, 



= -V • to', [J7] = if z = 0. 



Since to' is a Lipschitz vector field in fi (so, V • to' G L°°(K 2 )), by standard regularity theory 
for elliptic PDEs, we know that U is continuous in R 3 . We also may deduce some symmetry 
properties of U: First of all, the uniqueness of the stray field potential U G BC in (|93|) yields 
U(x, z) = U(x, —z) for every (x, z) el 3 . Also, remark that our vector field to' is anti-symmetric 
with respect to the origin, i.e., m'(x) — —m'(—x) which yields V • m'(x) = V • m'(—x) in R 2 . 
Again, by the uniqueness of the stray field potential U G BC, we deduce that U is symmetric in 
the in-plane variables with respect to the origin, i.e., 

U(x, z) = U(—x, ±z) for every (x, z) G R 3 . 

Also, the vector field to satisfies the symmetry relation m'{x) = (mi, —rri2)(—xi, X2) in R 2 , so that 
V • m'(x) — —V • m'(— X\,X2) in R 2 . It implies that 



U(xi,x 2 , z) = —U(—Xi,x 2 ,±z) for every (x, z) G 
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Similarly, U(xi,X2, z) = — U(x%, — X2, ±z) for every (x, z) € R 3 . In particular, it yields U(0, X2, z) = 
U(xi,0, z) = for every (x, z) e R 3 . 

In what follows, we compute upper bounds for / U(x, 0)V • m'(x)dx in several steps corre- 

Jn 

sponding to each subdomain of tt. In Qi t i U ( — U u\, m! is of vanishing divergence, therefore 

/ U{x 7 0)\7 ■m(x)dx ^0. (94) 

Jfil,lU(-fii,i)Uwi 

In the next step, we estimate 

[/(a:, 0)V ■m'(x)dx (95) 



where 

uj = n 1)2 U (-Oi, 2 ) U 2 ,i Uw 2 U w 3 . 
In the last step, we compute J B U(x, 0)V • m'(x) dx. 

Step 8: Upper bound for (|95p . The computation will be done according to the decomposition: 

V • m' = mi + m2 . In order to estimate / U(x, 0) m2 (x) dx, we use the following 

0x2 Jni, a u(-sii,z) dx 2 

argument (see also Proposition 3 in [8]): 



Lemma 3 Let L > 0, U : R 2 — > R and v : R — > R 6e suc/i £/ia£ v{x\) = v(—L) = v(L) for every 
\ x i I > L. Then 

( ^ 0) £ (Xl) " W) ( J R2 |VL/|2 dx ) ' 

where 

IMIii/a/m = ^ m in{ / |W| 2 dx : V(a?i,0) = v{x x ) for every x x € R}. (96) 
Here, we denote the homogeneous if ^-seminorm of v by 

:il^l 2 (£R, 



\/27r Ik 



where J-v £ S'(R) stands for the Fourier transform of v (as a tempered distribution), i.e., 

e'^vix^dxi, V£eR. 



One can also write 

Iff Hs)-v(t)\ 2 



Mlifl/2( 



2n L L \s-t\- 



dsdt (97) 



(see e.g., [8]). Another remark is that for even functions v (i.e., v(xi) = u(— xi)), the following 
estimate of ||u||jji/2(R) can be obtained via ([TO]) by considering the radial extension V of v in R 2 
(le.,V(x) = v(\x\)): 



h W)^2 7 r j w| r w dr - (98) 

Observe that is a general characterization of the H 1 / 2 — trace of if 1 — functions and it is valid 
in any dimension. 
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Proof of Lemma[3]. W.l.o.g., we can assume that v(xi) = v(—L) = v(L) = for every |xi| > L. 
Then Parseval's identity and the Cauchy-Schwarz inequality yield: 

L dv \ 2 ( f dv 

U(x 1 ,0) — (x 1 )dx 1 ) = U(x 1 ,0) — (x 1 )dx 1 



dxi J \J R ' dxi 



2 



< \^mu(;om,)\ 2 d^)( / i6ii^)(6)r^ 



f96l 1 / f 

- 2 IH| W)UJ W| d \ 

□ 

Writing each X\ € (1, 1 + <5) as xi = 1 + St with i S (0, 1), the X2— section in 0^2 passing 
through x\ is given by 

I X1 = {-^J\-5H\ y/1 - 6H*), 

so that O12 = U te (o i i){xi} x 4 r Since m 2 (xi, •) takes the same value at the boundary dI Xl for 
every t e (0, 1), we have by ([98) that: 

ll m 2(a;i,-)ll|i/2 (R) - f ^ |x 2 ||^^(a;i,x 2 )| 2 da; 2 < |^ |x 2 || J^(xi,x 2 )| 2 dx 2 = 0(1), 

where the upper bound 0(1) does not depend on X\. Therefore, Lemma [3] yields: 

U(x,0)^(x)dx Xl =± +tS 2 f s( [ U(x,0)^-(x)dx 2 )dt 
n li2 u(-fi li2 ) ox 2 Jo \J Imi dx 2 J 

<C^(^J \(v,-^jU(x,z)\ 2 dxdz^ f . (99) 
We apply the same argument to estimate / U(x, 0) ^ m2 (x) dx. By (I5TJ|) . we already know that 

IIM*ir)||^ (R) = fOTa11 Me(2£,l)- 

We deduce via Lemma [3] that: 

l/(»,Q)^(a?)d!C= / f / E^M)^ 2 -^)^ dan 

1/2 



9^2 i2<5<|x 1 |<l V7-1 9x 2 

1 



< 



^ lf^~' ^jU{x,z)\ 2 dx 2 dz S j ||rn 2 (xi,-)|lHi/2 (R) dxi 
" + 0(1) ) 1/2 (i l ( V '^) t/( "' Z)|2 ^" Z ) 1/2 - (1 ° 0) 



I log r?j 
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When estimating the same quantity in W3, a similar computation to (|99[) leads to 

1/2 



U(x,0)^-(x)dx < CVs( [ \V (xz) U{x,z)\ 2 dxdz\ 
ox 2 \Jk? J 



(101) 



In W2, a slightly different argument is used to estimate the quantity: 



U{x,0)^-{x)dx< ( j rs r.n, ',/./■ 



< 



\U(x,0)\ 4 dx 
<C8 5 ^\\U(;0)\\h^ 

1961 



1/4 



1/4 



dx 2 
dip 



4/3 \ 3/4 

(fa; 



9x 2 



4/3 X 3/4 

dx 



< C5 b/i (^J \(v,-^jU(x,z)\ 2 dxdz 



1/2 



(102) 



f din ^ 

It remains to estimate / U(x,0)— — (x) dx. In the region near the boundary, i.e., £l>n{\/l ~ 5 2 < 
1^2 1 < 1}, the same argument as in (|102l) yields: 



/ r 



U(x,0)^P±(x)dx< ( / |£/(:z;,0)| 4 efe 



1/4 



n{ v / T^P'<|x2|<i} 
1/2 



dip 



dx\ 



4/3 x 3/4 

dx 



J^\(v,^JU(x,z)\ 2 dxdz^ 



(103) 



For the interior region, i.e., Q n {ja^l < \/l — 5 2 }, we notice that e on f2 2 ,i and up to 



a translation, 4^ coincides on Oi^ and W3. Therefore, it is enough to estimate (by the above 
argument) the quantity 



U{x,Q)^p-(x)dx<[ / |[/(x,0)| 4 cfa 

Lo 3 n{\x 2 \<^T^} OX\ 



1/4 



< 



\U(x,0)\ 4 dx 



w 3 n{|a; 2 |<v / T I: 7 5 "} 

1/4 / r *JT=E I 



dip 



dx\ 



4/3 X 3/4 

dx 



-L( tta (a: 2 ) - A (x 2 )) 4 / 3 cfe 2 



The same computation as in (|83j) and (|84|) yields 



1 



^3 / at /3 (x 2 )dx 2 <5 2 /* 



1/5 



(t 2 + I) 2 / 3 



d< = 0(5 2 / 3 ) 



and 



^a/o ^ /3 (^)^ 2 = Q( J1/3|logA|4/3 ). 



Therefore, we deduce that 



U(x,0)^(x)dx<o( — 

I in{|2 2 |<v / I=^} OX! |fo; 



k^iUHhT- (im) 



Summing 0, (fTUO)) , (|TUT|) . ([TU2"]) . (fTU3")l and (fTUIl) , we obtain the following estimate for the stray 
field energy in J7: 

1/2 



Z7(x,0)V • to' cfe< 



7T + (1 )V /2 

I log?? I 



(105) 
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Step 9: Conclusion. It remains to estimate the stray field energy in B s as in (|102l) : 

/ s _ v { *. o)v . w ( / e> I (v, I) ci") " 2 ( J l Vm 'l 4/a ) 

^(IlKM" 2 - (106) 

By (|94|) , (| 105[) and (|106[) . we conclude that the total stray field energy is bounded by: 

7T + 0(1) 



l0gT/| 



i.e., by O, 



- / M^V-m') 
?7 Jr2 



2 , 2tt + o(1) , s 

dx < \> . (107) 

77|l0g7?| 

Finally, we estimate the last term of our energy given by the 7713— component. For our configuration 
ra, the only region in VL where m is not in-plane corresponds to the vortex core B e . There we have 



4r / m\dx^^r ( cos 2 (^-r) dx — 0(1). 



Combining with (|92|) and (1107j) . the conclusion follows. Remark that the constructed configuration 
m € H 1 ^, S 2 ) is only continuous. By the density of C 1 (ri, S 2 ) vector fields satisfying ([1]) in the 
space of i/ 1 (S7, 5 12 ) vector fields with ^ (for C 1 ' domains), one can smooth the configuration m 
so that the previous upper bound remains true. □ 

8 Appendix 

As mentioned in introduction, condition ([1} is necessary for a configuration to have finite stray-field 
energy in our model. To simplify the notation, we prove the statement for the case where d£l is a 
straight line: 

Proposition 2 Let f2 = (—00, 0)xl and m! G i/ 1 (f2,R 2 ). With the convention m! := m'ln, then 

2 

|V|" 1/2 (V • m!) dx < 00 implies that m^O, •) = in H 1/2 (R). 



Proof. We will show that 

m 1 (0,x 2 )(p(x 2 )dx 2 = for every <p e C™(R). (108) 

(Here, toi(0, •) represents the iJ 1 / 2 (R)— trace on the vertical line {x\ = 0} of mi € iJ 1 (J7,]R)). For 
a small e > 0, let £ e mimic the normal component of a Neel wall transition on a scale e with the 
size of the core of order e 2 : 



log j /, , 

if W\< J^7\ 



Ce(£l) = { 

elsewhere 
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We claim that (|108[) is equivalent to 

lim / mi(xi, X2)(p(x2) -r-^-(x\) dx\dx2 = for every (p £ C° 
e^0j R2 dxi 



(109) 



Indeed, we have: 



mi(0,x 2 )tp(x 2 ) dx^ 



< 



mi(xi,x 2 )ip(x 2 )- — (xi) dx\dx 2 - 
ax i 

mi(xi, X2)<p(x2)-^-(xi) dxidx 2 ~ I ^-(xi) I mi(0, x 2 )ip{x 2 ) dx 2 dx 



dxi 



dx\ 



L ^ {xi) i^KL ^ {s ' x2)ds ) dx2dxi 



< / \<p{x 2 )\ 



9mi 



dx i 



s,x 2 ) 



dsdx? 



IvK^III^ — {■,x 2 )\\ L 2^ ) dx 2 < s/s\\(p\\ L 2( 
ax i 



. dmi 
1 ~dx~[ 



\L 2 (n) 



(where we used that £ £ is increasing on R_ and J R j^(xi) dx\ = 1). In order to prove (|109p . we 
set ip(xi,x 2 ) = ( £ (xi)(p(x 2 ) and we write 

f f f dtp 

/ m\(x\,X2)tp{x2)—-^-{xi) dx\dx 2 = / m' ■ Vij) dx\dx 2 - I m 2 (xi,x 2 )Q(xi)- — (x 2 )dxidx 2 . 
Jr2 dxi J R 2 J M 2 dx 2 

Integrating by parts, we estimate the second term in the above RHS: 



m 2 (x 1 ,x 2 )( £ (xi)-^-(x 2 )dx 1 dx 2 
dx 2 



Ce(xi) [ ip(x 2 )^^-(x 1 ,x 2 )dx 2 dx 1 
Jm ox 2 

< IMIz, 2 (R) / Ce(^i)ll^^(a:i,OIU 2 (R)^i 



< Ve|MU 2 ( 



dm 2 
dx 2 



i 2 (n) 



(since Z e ( 2 e {x 1 )dx 1 "<~ n& /o^^log^s 2 +e 2 )d.s = 0( n ^)). The first term in the 
above RHS is estimated by interpolation: 



/ m' ■ V-0 




/ V • m'tp 


</ 


JR 2 




Jr 2 


7r2 



|Vr 1/2 (V-m') 



In order to conclude, we need to prove that \\ip\\ #1/2 (r2) — ► as e — >• 0. For that, we use (|96j) (valid 
in any dimension) for the following extension V : R 3 — > R of ^ given by V(x\, x 2 , z) = tp(r, x 2 ) = 
Q(r)if(x 2 ) for every (x\, x 2 , z) £ R 3 and r = \J x\ + z 2 : 



w 



I — I 

' dz ' 



dip 
dx 2 



and 



dx < 



2tt 



VV 



2 ,av, a 

' dz ' 



V? 2 (x 2 ) 



dxdz 



dr 



(r) 



-|*|J, mJ C^(r)* + M 

2 11 dl P 1,2 ' 



l80l 

< CI £ 



I logs] 



IMI 



— ^-(r) <ir 
dr 

as e ->• 0. 
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